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Abstract 

Guided by a geometric understanding developed in earlier works of Wang and Young, we 
carry out some numerical studies of shear-induced chaos. The settings considered include 
periodic kicking of limit cycles, random kicks at Poisson times, and continuous-time driving 
by white noise. The forcing of a quasi-periodic model describing two coupled oscillators is 
also investigated. In all cases, positive Lyapunov exponents are found in suitable parameter 
ranges when the forcing is suitably directed. 

Introduction 

This paper presents a series of numerical studies which investigate the use of shear in the pro- 
duction of chaos. The phenomenon in question can be described roughly as follows: An external 
force is applied to a system with tame, nonchaotic dynamics. If the forcing is strategically applied 
to interact with the shearing in the underlying dynamics, it can sometimes lead to the folding of 
phase space, which can in turn lead to positive Lyapunov exponents for a large set of initial condi- 
tions. This phenomenon, which we call shear-induced chaos, occurs in a wide variety of settings, 
including periodically-forced oscillators. For a topic as general as this, it is difficult to compile a 
reasonable set of references. We have not attempted to do that, but mention that the first known 
observation of a form of this phenomenon was by van der Pol and van der Mark 80 years ago [21]. 
Other references related to our work will be mentioned as we go along. 

The starting point of the present work is a series of papers by Wang and Young [|22ll23l[24ll25l . 
In these papers, the authors devised a method for proving the existence of strange attractors and 
applied their techniques to some natural settings; two of their main examples are periodically- 
kicked oscillators and systems undergoing Hopf bifurcations. They identified a simple geometric 
mechanism to explain how the chaotic behavior comes about. Because of the perturbative nature 
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of their analysis, however, the kicks in these results have to be followed by very long periods of 
relaxation. In other words, the chaos in these attractors develop on a very slow time scale. Relevant 
parts of the results of Wang and Young are reviewed in Sect. [TJ 

The aim of the present paper is to study shear-induced chaos in situations not accessible by 
current analytic tools. We believe that this phenomenon is widespread, meaning it occurs for large 
sets of parameters, and that it is robust, meaning it does not depend sensitively on the type of 
forcing or even background dynamics as long as certain geometric conditions are met. We validate 
these ideas through a series of numerical studies in which suitable parameters are systematically 
identified following ideas from ll23l and ll24ll . Four separate studies are described in Sects. [21- El 
The first three studies involve an oscillator driven by different types of forcing (both deterministic 
and stochastic); in these studies, the unforced system is a simple linear shear flow model. In the 
fourth study, the unforced dynamics are that of a coupled oscillator system described by a (periodic 
or quasi-periodic) flow on the 2-torus. 

The linear shear flow used in Studies 1-3 has been studied independently in ETl and [|23l . 
It is the simplest system known to us that captures all the essential features of typical oscillator 
models relevant to shear-induced chaos. Moreover, these features appear in the system in a way 
that is easy to control, and the effects of varying each are easy to separate. This facilitates the 
interpretation of our theoretical findings in more general settings in spite of the fact that numerical 
studies necessarily involve specific models. 

We mention that our results on shear flows are potentially applicable to a setting not discussed 
here, namely that of the advection and mixing of passive scalar tracers in (weakly compressible) 
flows. 

Finally, we remark that this work exploits the interplay between deterministic and stochastic 
dynamics in the following way: The geometry in deterministic models are generally more clear- 
cut. It enables us to extract more readily the relationship between quantities and to deduce the type 
of results these relationships may lead. Results for stochastic models, on the other hand, tend to be 
more provable than their counterparts in deterministic models, where competing scenarios lead to 
very delicate dependences on parameters. Our numerical results on stochastic forcing in Studies 
2-4 point clearly to the possibility of (rigorous) theorems, some versions of which, we hope, will 
be proved in the not too distant future. 

1 Rigorous Results and Geometric Mechanism 

In this section, we review some rigorous results of Wang and Young (mainly [|23ll24l . also Il22ll25ll ) 
and the geometric mechanism for producing chaos identified in the first two of these papers. We 
will focus on the case of limit cycles, leaving the slightly more delicate case of supercritical Hopf 
bifurcations to the reader. The material summarized in this section form the starting point for the 
numerical investigations in the present paper. 
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1.1 Strange Attractors from Periodically-Kicked Limit Cycles 

Consider a smooth flow $f on a finite dimensional Riemannian manifold M (which can be W^), and 
let 7 be a hyperbolic limit cycle, i.e. 7 is a periodic orbit of $t with the property that if we linearize 
the flow along 7, all of the eigenvalues associated with directions transverse to 7 have strictly 
negative real parts. The basin of attraction of 7, -6(7), is the set {x G M : ^ 7 as t ^ 00}. 

It is well known that hyperbolic limit cycles are robust, meaning small perturbations of the flow 
will not change its dynamical picture qualitatively. 

A periodically-kicked oscillator is a system in which "kicks" are applied at periodic time in- 
tervals to a flow $i with a hyperbolic limit cycle. For now let us think of a "kick" as a mapping 
K : M — M. If kicks are applied T units of time apart, then the time evolution of the kicked 
system can be captured by iterating its time-T map = $r o k. If there is a neighborhood U of 
7 such that k{U) C -6(7), and the relaxation time is long enough that points in k{U) return to U, 
i.e., FriU) C U, then F = nn>oF^{U) is an attractor for the periodically kicked system Ft. In a 
sense, F = F(k, T) is what becomes of the limit cycle 7 when the oscillator is periodically kicked. 
Since hyperbolic limit cycles are robust, F is a slightly perturbed copy of 7 if the kicks are weak. 
We call it an "invariant circle." Stronger kicks may "break" the invariant circle, leading to a more 
complicated invariant set. Of interest in this paper is when F is a strange attractor, i.e., when the 
dynamics in U exhibit sustained, observable chaos. 

Two theorems are stated below. Theorem 1 is an abstract result, the purpose of which is to 
emphasize the generality of the phenomenon. Theorem 2 discusses a concrete situation intended 
to make transparent the relevance of certain quantities. Let Leb(-) denote the Lebesgue measure 
of a set. 

Theorem 1. fl% Let $t be a flow with a hyperbolic limit cycle 7. Then there is an open set of 
kick maps K, with the following properties: For each k G /C, there is a set IS. = A{k) C M"*^ with 
Leb(A) > such that for each T ^ A, T is a "strange attractor" of Ft. 

The term "strange attractor" in the theorem has a well-defined mathematical meaning, which 
we will discuss shortly. But first let us take note of the fact that this result applies to all systems 
with hyperbolic limit cycles, independent of dimension or other specifics of the defining equations. 
Second, we remark that the kicks in this theorem are very infrequent, i.e. T ^ 1, and that beyond 
a certain Tq, the set A is roughly periodic with the same period as the cycle 7. 

The term "strange attractor" in Theorem 1 is used as short-hand for an attractor with a package 
of well defined dynamical properties. These properties were established for a class of rank-one 
attractors (see ll22l for the 2-dimensional case; a preprint for the n-dimensional case will appear 
shortly). In |[22l |25l , the authors identified a set of conditions that implies the existence of such 
attractors, and the verification of the conditions in [|22l |25l in the context of Theorem 1 is carried 
out in |[24]| (see also [fTTl and IfTTl for other applications of these ideas). We refer the reader to the 
cited papers for more details, and mention only the following three characteristics implied by the 
term "strange attractor" in this section. 

(1) There is a set V of full Lebesgue measure in the basin of attraction of F such that orbits 
starting from every x G V have (strictly) positive Lyapunov exponents. 
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Figure 1 : The stretch-and-fold action of a kick followed by relaxation in the presence of shear. 
(2) Ft has an ergodic SRB measure ji, and for every continuous observable y?, 



if dfx as 77. — >• oo for every X G V. 



(3) The system (FT,fi) is mixing; in fact, it has exponential decay of correlations for Holder 
continuous observables. 

An important remark before leaving Theorem 1 : Notice that the existence of "strange attrac- 
tors" is asserted for Ft for only a positive measure set of T, not for all sufficiently large T. This is 
more a reflection of reality than a weakness of the result. For large enough T in the complement 
of A, the attractor is guaranteed to contain horseshoes, the presence of which will lead to some 
semblance of chaotic behavior. After a transient, however, typical orbits may (or may not) tend 
to a stable equilibrium. If they do, we say Ft has transient chaos. This is to be contrasted with 
properties (l)-(3) above, which represent a much stronger form of chaos. 

The next result has an obvious analog in n-dimensions (see (241), but the 2-D version illustrates 
the point. 

Theorem 2. |!24| Consider the system 

e = 1 + ay 

y = .Xy + A-H{e).Y.n=,5{t-nT) 

where [6, y) G 5^ x R are coordinates in the phase space, X,a,A > are constants, and H : 
^ is a nonconstant smooth function. If the quantity 

^ A shear , . , , „ 

— ■ A = ■ kick amplitude 

A contraction rate 

is sufficiently large (how large depends on the forcing function H), then there is a positive measure 
set A C M"*" such that for all T G A, has a strange attractor in the sense above. 

Here, the term involving H(9) defines the kick, and 7 = x {0}. We explain intuitively the 
significance of the quantity jA. As noted earlier, to create a strange attractor, it is necessary to 
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Figure 2: Geometry of folding in relation to the VF^* -foliation. Shown are the kicked image of a segment 
7o and two of its subsequent images under 

"break" the limit cycle. The more strongly attractive 7 is, the harder it is to break. From this we 
see the advantage of having A small. By the same token, a stronger forcing, i.e., larger A, helps. 
The role of a, the shear, is explained pictorially in Fig. [TJ Since the function H is required to be 
nonconstant, let us assume the kick drives some points on the limit cycle 7 up and some down, as 
shown. The fact that a is positive means that points with larger y-coordinates move faster in the 9- 
direction. During the relaxation period, the "bumps" created by the kick are stretched as depicted. 
At the same time, the curve is attracted back to the limit cycle. Thus, the combination of kicks 
and relaxation provides a natural mechanism for repeated stretching and folding of the limit cycle. 
Observe that the larger the differential in speed in the 6'-direction, i.e. the larger a, and the slower 
the return to 7, i.e. the smaller A, the more favorable the conditions are for this stretch-and-fold 
mechanism. 

1.2 Geometry and Singular Limits 

In Eq. ([U), the quantities A, a and A appear naturally. But what about in general limit cycles, where 
the direction of the kicks vary? What, for example, will play the role of a, or what we called shear 
in Eq. ([T])? The aim of this subsection is to shed light on the general geometric picture, and to 
explain how the dynamics of for large T can be understood. 

Geometry of Ft and the Strong Stable Foliation 

Let 7 be a hyperbolic limit cycle as in the beginning of Sect. 1.1. Through each x G 7 passes 
the strong stable manifold of x, denoted W^'^{x) IfTOl . By definition, W^'^{x) = {y E M : 
(i($t(|/), $t(x)) — i> as t ^ 00}; the distance between and ^tiv) in fact decreases ex- 

ponentially. Some basic properties of strong stable manifolds are: (i) W^^^x) is a codimension 
one submanifold transversal to 7 and meets 7 at exactly one point, namely x; (ii) $j(iy*(a;)) = 
W'^'^{^t{x)), and in particular, if the period of 7 is p, then ^p{W'^'^{x)) = W^^{x); and (iii) the 
collection {W^^{x),x G 7} foliates the basin of attraction of 7, that is to say, they partition the 
basin into hypersurfaces. 

We examine next the action of the kick map k, in relation to iy*''-manifolds. Fig. [2lis analogous 
to Fig. [11 it shows the image of a segment 70 of 7 under Ft = $t o For illustration purposes, we 
assume 70 is kicked upward with its end points held fixed, and assume T = np for some n G Z+ 
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(otherwise the picture is shifted to another part of 7 but is qualitatively similar). Since $„p leaves 
each "^-manifold invariant, we may imagine that during relaxation, the flow "slides" each point 
of the curve ^(70) back toward 7 along W^^*-leaves. In the situation depicted, the effect of the 
folding is evident. 

Fig. I2] gives considerable insight into what types of kicks are conducive to the formation of 
strange attractors. Kicks along iy**-leaves or in directions roughly parallel to the iy**-leaves will 
not produce strange attractors, nor will kicks that essentially carry one H^**-leaf to another. What 
causes the stretching and folding is the variation in how far points x G 7 are moved by k as 
measured in the direction transverse to the 14^ '"'-leaves. Without attempting to give a more precise 
characterization, we will refer to the type of chaos that results from the geometry above as shear- 
induced chaos. We emphasize that the occurrence of shear-induced chaos relies on the interplay 
between the geometries of the kicks and the underlying dynamical structures. 

Returning to the concrete situation of Theorem 2, since Eq. ([T]) without the kick term is linear, 
it is easy to compute strong stable manifolds. In {9, i/)-coordinates, they are lines with slope — A/cr. 
Variations in kick distances here are guaranteed by the fact that H is nonconstant. With H fixed, it 
is clear that the larger a/X and A, the greater these variations. Note that the use of the word kick 
"amplitude" in the statement of Theorem 2 is a little misleading, for it is not the amplitude of the 
kicks per se that leads to the production of chaos. 

Singular Limits of as T — i> 00 

When T 1, i.e. when kicks are very infrequent, the map Ft sends a small tube Ut around 7 
back into itself. This is an example of what is called a rank-one map in [25 1. Roughly speaking, 
a rank-one map is a smooth map whose derivative at each point is strongly contractive in all but 
one of the directions. Rank-one maps can be analyzed using perturbative methods if they have 
well-defined "singular limits." In the context of limit cycles, these singular limits do exist; they are 
a one-parameter family of maps {fa : 7 0} obtained by letting T ^ 00 in the following way: For 
each a E [0,p) (recall that p = period of 7), let 

fa{x) := lim for aU a: G 7 . (2) 

n— >oo 

Equivalently, fa{x) is the unique point y E j such that k{x) E W'^^{y). Notice that faix) = 
fo{x) + a (mod 1), where we identify 7 with [0, 1] (with the end points identified). For Eq. O, /„ 
is easily computed to be 

fai9) = 9 + a + jA-H{9), (3) 

where the right side should again be interpreted as mod 1. (In the setting of driven oscillators, 
singular limits are sometimes known as "phase resetting curves"; they have found widespread use 
in e.g. mathematical biology [|26l l9ll.) 

It is shown in ll22l l23l l24l l25l that a great deal of information on the attractor T of Ft for 
T 2> 1 can be recovered from these singular limit maps. The results are summarized below. These 
results hold generally, but as we step through the 3 cases below, it is instructive to keep in mind 
Eq. ([T]) and its singular limit ([3]), with jA increasing as we go along: 
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(i) If fa is injective, i. e. , it is a circle diffeomorphism, the attractor T for Ft is an invariant circle. 
This happens when the kicks are aimed in directions that are "unproductive" (see above), or 
when their effects are damped out quickly. In this case, the competing scenarios on T are 
quasi-periodicity and "sinks," i.e. the largest Lyapunov exponent of Ft is zero or negative. 

(ii) When fa loses its injectivity, the invariant circle is "broken". When that first happens, the 
expansion of the 1-D map fa is weak, and all but a finite number of trajectories tend to sinks. 
This translates into a gradient type dynamics for Ft. 

(iii) If fa is sufficiently expanding away from its critical points, T contains horseshoes for all 
large T. For an open set of these T, the chaos is transient, while on a positive measure 
set. Ft has a strange attractor with the properties described in Sect. 1.1. These are the two 
known competing scenarios. (They may not account for all T.) Since Ft ~ Ft+up for large 
T, both sets of parameters are roughly periodic. 

The analyses in the works cited suggest that when horseshoes are first formed, the set of 
parameters with transient chaos is more dominant. The stronger the expansion of /„, the 
larger the set of parameters with strange attractors. In the first case, the largest Lyapunov 
exponent of Ft may appear positive for some time (which can be arbitrarily long) before 
turning negative. In the second case, it stays positive indefinitely. 

1.3 Limitations of Current Analytic Techniques 

In hyperbolic theory, there is, at the present time, a very large discrepancy between what is thought 
to be true and what can be proved. Maps that are dominated by stretch-and-fold behavior are 
generally thought to have positive Lyapunov exponents - although this reasoning is also known 
to come with the following caveat: Maps whose derivatives expand in certain directions tend to 
contract in other directions, and unless the expanding and contracting directions are well separated 
(such as in Anosov systems), the contractive directions can conspire to form sinks. This is how 
the transient chaos described in Sect. 1 .2 comes about. Still, if the expansion is sufficiently strong, 
one would expect that positive Lyapunov exponents are more likely to prevail - even though for 
any one map the outcome can go either way. Proving results of this type is a different matter. Few 
rigorous results exist for systems for which one has no a priori knowledge of invariant cones, and 
invariant cones are unlikely in shear-induced chaos. 

The rigorous results reviewed in the last two subsections have the following limitations: (i) 
They pertain to Ft for only very large T. This is because the authors use a perturbative theory 
that leans heavily on the theory of 1-D maps. No non-perturbative analytic tools are currently 
available, (ii) A larger than necessary amount of expansion is required of the singular limit maps 
fa in the proof of strange attractors. This has to do with the difficulty in locating suitable parameters 
called Misiurewicz points from which to perturb. (This problem can be taken care of, however, by 
introducing more parameters.) We point out that (i) and (ii) together exacerbate the problem: fa 
is more expanding when A is small, but if Ft = $t o «^ is to be near its singular limit, then e~^^ 
must be very small, i.e. XT must be very large. 
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That brings us to the present paper, the purpose of which is to supply numerical evidence to 
support some of our conjectured ideas regarding situations beyond the reach of the rigorous work 
reviewed. Our ideas are based on the geometry outlined in Sect. 11.21 but are not limited to periodic 
kicks or to the folding of limit cycles. 



2 Study 1: Periodically-Kicked Oscillators 

Our first model is the periodic kicking of a linear shear flow with a hyperbolic limit cycle. The 
setting is as in Theorem 2 with H{9) = sin(27r^), i.e., we consider 

e = 1 + ay, 

y = -Xy + A- sin(27r^) ■ ^^=0 ^(^ ~ ^T) , ^ ^ 

where (9, y) E x M, = [0, 1] with the two end points of [0, 1] identified. In the absence of 
kicks, i.e., when A = 0, ^t{z) tends to the limit cycle 7 = S*^ x {0} for all 2 e x R. As before, 
the attractor in the kicked system is denoted by T. The parameters of interest are: 

a = amount of shear, 

A = damping or rate of contraction to x {0}, 
A = amplitude of kicks, and 
T = time interval between kicks. 

Our aim here is to demonstrate that the set of parameters with chaotic behavior is considerably 
larger than what is guaranteed by the rigorous results reviewed in Sect. [H and to gain some insight 
into this parameter set. By "chaotic behavior," we refer in this section to the property that Ft has 
a positive Lyapunov exponent for orbits starting from a "large" set of initial conditions, i.e. a set 
of full or nearly full Lebesgue measure in the basin of attraction of T. More precisely, we assume 
that such Lyapunov exponents are well defined, and proceed to compute the largest one, which we 

call Ajnax- 

We begin with some considerations relevant to the search for parameters with Amax > 0: 

(a) It is prudent, in general, to ensure that orbits do not stray too far from 7. This is because 
while the basin of attraction of 7 in this model is the entire phase space, the basin is bounded 
in many other situations. We therefore try to keep T C {\y\ < b} with relatively small b. 
This is guaranteed if A is small enough that e~^^{b + A) < b; the bound is improved if, for 
example, no point gets kicked to maximum amplitude two consecutive iterates. 



(b) Let {6t, yr) = Ft{6o, y^). A simple computation gives 

I - P-^'^\ rmnH n 

(5) 



Ot = 9o + T + ^-[yo + A sin(27r^o)] ■ (1 - e"^^) (mod 1) 
Vt = e-^'^[yo + Asm{2TTeo)] . 



For b relatively small, we expect the number ^(1 — e~^^) to be a good indicator of chaotic 
behavior: if it is large enough, then Ft folds the annulus {\y\ < b} with two turns and maps 
it into itself. The larger this number, the larger the folds, meaning the more each of the 
monotonic parts of the image wraps around in the ^-direction. 



8 



Summary of Findings. 

( i) With the choice of parameters guided by (a) and (b) above, we find that as soon as the folding 
described in (b) is definite, Ft becomes "possibly chaotic", meaning A^ax is seen numer- 
ically to oscillate (wildly) between positive and negative values as T varies. We interpret 
this to be due to competition between transient and sustained chaos; see (Hi) in Sect. 1.2. 
For larger jA, i.e., as the stretching is stronger, and for T beyond an initial range, this 
oscillation stops and Amax becomes definitively positive for all the values ofT computed. 

(ii) As for the range of parameters with chaotic dynamics, we find that A^ax > occurs under 
fairly modest conditions, e.g., for jA = 3, we find Amax > starting from about T ^ 3, 
which is very far from the "T oo" in rigorous proofs. Also, while shear-induced chaos is 
often associated with weak damping, we find that the phenomenon occurs as well for larger 
X, e.g., for A ~ 1, provided its relation to the other parameters are favorable. 

Supporting Numerical Evidence. Figures |3] and |4] show the largest Lyapunov exponent Amax of 
Ft versus the kick period T. (Note that this is the expansion rate per kick period and is T times 
the rate per unit time.) In Fig.[3l A and A are fixed, and a is increased. We purposefully start with 
too small a a so that we may see clearly the gradual changes in A^ax- The results are in excellent 
agreement with the description at the end of Sect. 11.21 (which pertains to regimes with very large 
T), even though T is not so large here: In the top picture, where jA is small, the plot confirms a 
competition between quasi-periodicity and sinks; in the middle picture, we see first A^ax becoming 
increasingly negative, then transitions into a competition between transient and sustained chaos, 
with the latter dominating in the bottom picture. Fig.|4]shows the same phenomena in reverse order, 
with a and A fixed and A increasing. Notice that even for a, A and A leading to chaotic dynamics, 
Amax is negative for small T. This is in agreement with the influence of the factor (1 — e~'^^) in 
Eq.©. 

As explained in (a) above, when AT is too small relative to A, orbits stray farther from 7. 
Data points corresponding to parameters for which this happens are marked by open squares. For 
purposes of demonstrating the phenomena in question, there is nothing wrong with these data 
points, but as explained earlier, caution must be exercised with these parameters in systems where 
the basin of 7 is smaller. 

Simulation Details. The numbers A^ax are computed by iterating the map in Eq. ([5]) and its Jaco- 
bian, and tracking the rate of growth of a tangent vector. We use 4 x 10^ iterates of Ft in each 
run. Mindful of the delicate situation due to competition between transient and sustained chaos, 
and to lower the possibility of atypical initial conditions, we perform 10 runs for each choice of 
(cr, A, A, T), using for each run an independent, random (with uniform distribution) initial condi- 
tion (6^0, yo) G [0, 1) X [—0.1, 0.1]. Among the 10 values of Amax computed, we discard the largest 
and the smallest, and plot the maximum and minimum of the remaining 8. As one can see in Figs.|3] 
andlH the two estimates occasionally do not agree. This may be because not all initial conditions 
in the system have identical Lyapunov exponents, or it may be that the convergence to the true 
value of Amax is sufficiently slow and more iterates are needed, i.e. there are long transients. 
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Figure 3: Effect of increasing shear on the Lyapunov exponents of the periodically-kicked linear shear flow. 
Squares indicate that the corresponding F^^-orbit has veered outside the region \y\ < 0.15. Upper and lower 
estimates of Amax are both shown (see Simulation Details). 
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Figure 4: Effect of increasing damping on the Lyapunov exponents of the periodically-kicked linear shear 
flow. Squares indicate that the corresponding FT-orhit has veered outside the region \y\ < 0.15. 
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3 Study 2: Poisson Kicks 



We consider next a variant of Eq. dH) in which deterministic, periodic kicks are replaced by "ran- 
dom kicks." Here, random kicks refer to kicks at random times and with random amplitudes. More 
precisely, we consider 

e=l + ay (6) 
y = -\y + sin(27r^) ^ AJ{t - T„) 

n 

where the kick times T„ are such that T,i+i — T„, n = 0, 1, 2, ■ ■ ■ , are independent exponential 
random variables with mean T, and the kick amplitudes A„ are independent and uniformly dis- 
tributed over the interval [0.8 A, 1.2 A] for some A > 0. (We do not believe detailed properties 
of the laws of T and A have a significant impact on the phenomena being addressed.) The analog 
here of the time-T map in Study 1 is the random map F = $t ° Ka where T and A are random 
variables. 

By the standard theory of random maps, Lyapunov exponents with respect to stationary mea- 
sures are well defined and are nonrandom, i.e. they do not depend on the sample path taken [12J. 
Notice that if a 7^ 0, the system ^ has a unique stationary measure which is absolutely continu- 
ous with respect to Lebesgue measure on x R: starting from almost every zq E x M, after 
one kick, the distribution acquires a density in the y-direction; since vertical lines become slanted 
under due to cr 7^ 0, after a second kick the distribution acquires a (two-dimensional) density. 

In terms of overall trends, our assessment of the likelihood of chaotic behavior follows the 
analysis in Study 1 and will not be repeated. We identify the following two important differences: 

(a) Smooth dependence on parameters. Due to the averaging effects of randomness, we expect 
Lyapunov exponents to vary smoothly with parameter, without the wild oscillations in the 
deterministic case. 

(b) Effects of large deviations. A large number of kicks occurring in quick succession may have 
the following effects: 

(i) They can cause some orbits to stray far away from 7 = S*^ x {0}. This is guaranteed 
to happen, though infrequently, in the long run. Thus, it is reasonable to require only 
that a large fraction — not all — of the stationary measure (or perhaps of the random 
attractors T^) to lie in a prescribed neighborhood of 7. 

(ii) It appears possible, in principle, for a rapid burst of kicks to lead to chaotic behavior 
even in situations where the shear is mild and kick amplitudes are small. To picture this, 
imagine a sequence of kicks sending (or maintaining) a segment far from 7, allowing 
the shear to act on it for an uncharacteristically long time. One can also think of such 
bursts as effectively setting A to near temporarily, creating a very large ^A. On the 
other hand, if a is small, then other forces in the system may try to coax the system to 
form sinks between these infrequent events. We do not have the means to assess which 
scenario will prevail. 
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Figure 5: Lyapunov exponents for the linear shear flow with Poisson kicks. Squares indicate the correspond- 
ing orbit spends more than 20% of the time in the region \y\ > 0.1. 

Summary of Findings. In terms of overall trends, the results are consistent with those in Study 1. 
Two differences are observed. One is the rapid convergence o/ Amax and their smooth dependence 
on parameters. The other is that positive Lyapunov exponents for F are found both for smaller 
values of jA and for apparently very small T (which is impossible for periodic kicks), lending 
credence to the scenario described in (b)(ii) above. 

Supporting Numerical Evidence. Fig.|5]shows Amax as a function of the mean kick interval T. As in 
Study 1, we first show the effects of increasing a and then the effects of increasing A. Without the 
oscillations seen previously, the present plots are straightforward to interpret. In case one wonders 
how Ajnax curves can switch from strictly-decreasing to strictly-increasing behavior, the middle 
panel of Fig. (S^b) catches such a swtich "in the act." Squares indicate that the orbit computed 
spends > 20% of its time outside of the region {\y\ < 0.1}. 
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4 Study 3: Continuous-Time Stochastic Forcing 



In this section, we investigate the effect of forcing by white noise. The resuhing systems are 
described by stochastic differential equations (SDEs). We consider two ways to force the system: 

Study 3a: Degenerate white noise applied in chosen direction: 

d9={l + ay) dt (7) 
dy = —Xy dt + a sin(27r6') dBt 

Study 3b: Isotropic white noise: 

de={l + ay) dt + a sin(27r^) dB} (8) 
dy = —\ydt + a sin(27r6') dB^ 

In Study 3a, Bt is standard 1-dimensional Brownian motion (meaning with variance = 1). In Study 
3b, {Bl, B^) is a standard 2-D Brownian motion, i.e., they are independent standard 1-D Brownian 
motions. For definiteness, we assume the stochastic terms are of Ito type. Notice that the two 
parameters A and T in Studies 1 and 2 have been combined into one, namely a, the coefficient of 
the Brownian noise. 

By standard theory [HI [IS, the solution process of an SDE can be represented as a stochastic 
flow of diffemorphisms. More precisely, if the coefficients of the SDE are time-independent, then 
for any time step At > 0, the solution may be realized, sample path by sample path, as the 
composition of random diffeomorphisms ■ ■ ■ 0/30/2°/!, where the /j are chosen i.i.d. with 
a law determined by the system (the /i are time-At flow-maps following this sample path). This 
representation enables us to treat an SDE as a random dynamical system and to use its Lyapunov 
exponents as an indicator of chaotic behavior. It is clear that system ^ has a unique invariant 
density, which is the solution of the Fokker-Planck equation. Even though the stochastic term in 
system (|7]) is degenerate, for the same reasons discussed in Study 3, it too has a unique stationary 
measure, and this measure has a density. The Lyapunov exponents considered in this section are 
with respect to these stationary measures. 

Before proceeding to an investigation of the two systems above, we first comment on the case 
of purely additive noise, i.e. Eq. ([8]) without the sin(27r6') factor in either Brownian term. In this 
case it is easy to see that all Lyapunov exponents are < 0, for the random maps are approximately 
time-At maps of the unforced flow composed with random (rigid) translations. Such a system is 
clearly not chaotic. 

With regard to system (|7]), we believe that even though the quantitative estimates from Study 
1 no longer apply, a good part of the qualitative reasoning behind the arguments continues to be 
valid. In particular, we conjecture that 

(a) trends, including qualitative dependences on a and A, are as in the previous two studies; 

(b) the effects of large deviations noted for Poisson kicks (Study 2, item (b)) are even more 
prominent here, given that the forcing now occurs continuously in time. 
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As for system ([8]), we expect it to be less effective in producing chaos, i.e. more inclined to form 
sinks, than system ([7]). This expectation is based on the following reasoning: Suppose first that we 
force only in the ^^-direction, i.e., suppose the dB'f term in ([8]) is absent. Then the stochastic flow 
leaves invariant the circle x {0}, which is the limit cycle of the deterministic part of the system. 
A general theorem tells us that when a random dynamical system on a circle has an invariant 
density, its Lyapunov exponent is always < 0; in this case, it is in fact strictly negative because 
of the inhomogeneity caused by the sine function lfT2]| . Thus the corresponding 2-D system has 
"random sinks." Now let us put the y-component of the forcing back into the system. We have 
seen from previous studies that forcing the y direction alone may lead to chaotic behavior. The 
tendency to form sinks due to forcing in the 6'-direction persists, however, and weakens the effect 
of the shear- induced stretching. 

We now discuss the results of simulations performed to validate these ideas. 

Summary of Findings. 

(i) In the case o/ degenerate white noise, the qualitative dependence of K^.^^ on a and A are as 
expected, and the ejfects of large deviations are evident. In particular, Amax is positive for 
very small values of a, X and a provided j is large. This cannot happen for periodic kicks; 
we attribute it to the effect of large deviations. 

( ii) Isotropic white noise is considerably less effective in producing chaos than forcing in the 
y -direction only, meaning it produces a smaller (or more negative) Amax- 

( Hi) In both cases, we discover the following approximate scaling: Under the scaling transforma- 
tions \ ^ kX,a ^ ka and a ^ \/ka, A^ax transforms approximately as Amax ^ ^Amax- 
In the case of degenerate white noise, when both a and - are not too small (e.g., > 3), this 
scaling gives excellent predictions of K^g^^for the values computed. 

We remark that (iii) does not follow by scaling time in the SDE. Indeed, scaling time by k in 
Eq. ([7]), we obtain 

de = {k + kay) dt , (9) 
dy = —kXy dt + v^asin(27r^) dBt . 

Thus the approximate scaling in (iii) asserts that the Lyapunov exponent of system Q, equivalently 
k times the Amax for Eq. (|7]), is roughly equal to that of the system obtained by changing the first 
equation in © to (i6' = (1 + kay)dt. In other words, Amax seems only to depend minimally on the 
frequency of the limit cycle in the unforced system. 

Supporting Numerical Evidence. Plots of Amax as functions of a are shown in Figs. [6]-[9l 

In Fig.|6l the forcing is degenerate, and for fixed a, Amax decreases with increasing damping as 
expected. Notice that compared to the two previous studies, a somewhat larger damping is required 
to maintain a good fraction of the attractor near 7. 
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Figure 6: Lyapunov exponents for the linear shear flow driven by degenerate white noise (Eq. Q). Open 
squares indicate that the corresponding orbits spend more than 20% of the time in the region \y\ > 0.3; 
shaded squares do the same for the region \y\ > 0.5. 
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Figure 7: Lyapunov exponents for the linear- sheai^ flow driven by degnerate white noise, for small values of 
a and A. 

Fig.|7]shows that A^ax is positive for values of a and A as small as 0.2 and 0.01, and white noise 
amplitudes a close to 0. Notice first that this is consistent with the scaling conjectured in (iii) above, 
and second that in the case of periodic kicks, comparable values of a and A would require a fairly 
substantial kick, not to mention long relaxation periods, before chaotic behavior can be produced. 
We regard this as convincing evidence of the significant effects of large deviations in continuous- 
time forcing. (It must be pointed out, however, that in our system, the basin of attraction is the 
entire phase space, and a great deal of stretching is created when |y| is large. That means system 
([7]) takes greater advantage of large deviations than can be expected ordinarily. 

Fig.[8ta) shows A^ax in the isotropic case for the same parameters as in Fig.[6l A comparison 
of the two sets of results confirms the conjectured tendency toward negative exponents when the 
forcing is isotropic. Fig.[8]^b) shows that this tendency can be overcome by increasing a. 

Fig. |9] shows four sets of results, overlaid on one another, demonstrating the scaling discussed 
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(b) a = 7.0 

Figure 8: Lyapunov exponent for the linear shear flow driven by isotropic white noise (Eq. Squares 
indicate that the corresponding orbits spend more than 20% of the time in the region \y\ > 0.3. 

in item (ill) above. Fixing ^ = 6, we show the graphs of Amax/o" as functions of a/^/a for four 
values of a. The top two curves (corresponding to cr = 6 and 9) coincide nearly perfectly. Similar 
approximate scalings, less exact, are observed for smaller values of ^, both when A^ax is positive 
and negative. 

Simulation Details. We compute Lyapunov exponents numerically by solving the corresponding 
variational equations (using an Euler solver with time steps of 10^^) and tracking the growth rate of 
a tangent vector. To account for the impact of the realization of the forcing on the computed expo- 
nents, for each choice of {a, A, a) we perform 12 runs in total, using 3 independent realizations of 
the forcing and, for each realization, 4 independent initial conditions (again uniformly-distributed 
in [0, 1) X [—0.1,0.1]). For almost all the parameter values, the estimates agree to fairly high 
accuracy, so we simply average over initial conditions and plot the result. 
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Figure 9: Evidence of scaling: We fix ^ and plot Amax/c as functions of the rescaled drive amplitude a/y^; 
from top to bottom, the curves are in order of decreasing a. 

Related Results. The asymptotic stability of dynamical systems driven by random forcing has 
been investigated by many authors using both numerical and analytic methods. Particularly rele- 
vant to our study are results pertaining to the random forcing of oscillators (such as Duffing-van 
der Pol oscillators) and stochastic Hopf bifurcations; see e.g. Il2l|3l|7l|5l[6l[8l[l8lll6l. Most of the 
existing results are perturbative, i.e., they treat regimes in which both the noise and the damping 
are very small. Positive Lyapunov exponents are found under certain conditions. We do not know 
at this point if the geometric ideas of this paper provide explanations for these results. 

5 Study 4: Sheared-Induced Chaos in Quasiperiodic Flows 

Model and Background Information 

In this section, we will show that external forcing can lead to shear-induced chaos in a coupled 
phase oscillator system of the form 

I{t) 



The governing equations are 

e^ = v^ + z{e^)[a^g{92) + I{t)], (10) 
62 = iy2 + z{e2)[aikg{ei)]. 

The state of the system is specified by two angles, (6*1, 6*2), so that the phase space is the torus = 
[0, 1)^. The constants ui and z/2 are the oscillators' intrinsic frequencies; we set z/i = 1 and 1^2 = 1.1 
(representing similar but not identical frequencies). The constants ag and Ofb govern the strengths 
of the feedforward and feedback couplings. The oscillators are pulse-coupled: the coupling is 
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mediated by a bump function g supported on [— ^, ^] and normalized so that g{6) dO = 1. 
The function z{0), which we take to be z{9) = 2^[1 — cos(27r6')], specifies the sensitivity of the 
oscillators to perturbations when in phase 6. Finally, we drive the system with an external forcing 
I{t), which is applied to only the first oscillator. This simple model arises from neuroscience 
Il26ll20ll and is examined in more detail in ifTSl . 

Let $t denote the flow of the unforced system, i.e., with I{t) = 0. Flowlines are roughly 
northeasterly and are linear except in the strips {|6'i| < ^} and {\92\ < ^}, where they are 
bent according to the prescribed values of afj and Ofb- Let p denote the rotation number of the 
first return map of $f to the cross-section {02 = 0}. It is shown in il4j that for cig = 1, p 
is monotonically increasing (constant on extremely short intervals) as one increases Ofb, until it 
reaches 1 at Ofb = ~ 1.4, after which it remains constant on a large interval. At Ofb = a^, a 
limit cycle emerges in which each oscillator completes one rotation per period; we say the system 
is 1:1 phase-locked, or simply phase-locked. In [[T4||. it is shown numerically that forcing the 
system by white noise after the onset of phase-locking leads to A^ax > 0. The authors of [13 
further cite Wang- Young theory (the material reviewed in Sect. [I]) as a geometric explanation for 
this phenomenon. 

In this section, we provide geometric and numerical evidence of shear-induced chaos both 
before and after the onset of phase-locking at Ofb = a^. Our results for Ofb > support the 
assertions in [14J. For Ofb < a^, they will show that limit cycles are not preconditions for shear- 
induced chaos. We will show that in Eq. (fTOl) . the mechanism for folding is already in place 
before the onset of phase-locking, where the system is quasi-periodic or has periodic orbits of very 
long periods; the distinction between these two situations is immaterial since we are concerned 
primarily with finite-time dynamics. In the rest of this section, we will, for simplicity, refer to the 
regime prior to the onset of phase-locking as "near-periodic." 

Folding: Geometric Evidence of Chaos 

The dynamical picture of kicks followed by a period of relaxation has a simpler, more clear-cut 
geometry than that of continuous, random forcing. Thus we use the former to demonstrate why 
one may expect chaotic behavior over the parameter ranges in question. The kick map is denoted 
by K as in Section 1 . 

Folding in the periodic (i.e. phase-locked) regime. We will use = 1.47 for illustration purposes; 
similar behavior is observed over a range of Ofb from 1.4 to 1.6. Note that the system is phase- 
locked for a considerably larger interval beyond afb = 1.6, but the strength of attraction grows with 
increasing Ofb, and when the attraction becomes too strong, it is harder for folding to occur. 

Fig. [To] shows the limit cycle 7 {blue curve) of the unforced system at Cfb = 1.47; more 
precisely, it shows a "lift" of 7 to M^, identifying the torus with M^/Z^. Also shown is the 
image ^(7) of the cycle after a single kick {red curve), where the kick map k corresponds to 
I{t) = AJ2n^i^ ~ ^'^) ^^^^ A = 1-5, i-e. K is given hy k, = lime^o Ke{e) where ^^(t) is the 
solution of ^1 = jz{6i),62 = 0. Notice the special form of the kicks: k acts horizontally, and does 
not move points on 9i = 0. In particular, n fixes a unique point (0, b) on the cycle; this point is, in 
fact, not affected by any kick of the form considered in Eq. (flOl) . Several segments of strong stable 
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Figure 10: The strong-stable foliation of the system (ITOl ) in the phase-locked regime. Blue: a lift of the limit 
cycle. Red: the image of the cycle after a single kick. Green: strong-stable foliation. Here, the parameters 
are ui = 1,V2 = 1.1, = 1, and Ofb = 1.47. 

manifolds {green curves) of the unforced system are drawn. Recall that if p is the period of cycle 
and n E Z+, then $„p(k(2;)) lies on the '"'-curve through k{z) and is pulled toward the cycle as 
n increases (see Sect. 11.21) . From the relation between the l^^'^-curves and the cycle, we see that 
for 2 G 7, $f (k(2;)) will lag behind ^tiz) during the relaxation period. Notice in particular that 
there are points on ^(7) above the line 62 = b that are pulled toward the part of 7 below 62 = b. 
Since (0, b) stays put, we deduce that some degree of folding will occur if the time interval between 
kicks is sufficiently long. 

Fig, nn illustrates how this folding happens through three snapshots. We begin with a segment 
7o C 7 between 62 = and 62 = 1 (blue curve) and its image after a single kick {red curve). Both 
curves are then evolved forward in time and their images at t = 2.5 and t = 3.5 are shown. The 
purple dot marks the point on 70 which does not move when kicked. Notice that these pictures are 
shown in a moving frame to emphasize the geometry of $4(^(70)) relative to $t(7o). 

Folding in the near-periodic regime. Fig. [12] shows snapshots of a similar kind for Ofb = 1.2; this 
value of Ofb puts the system in the near-periodic regime. The snapshots begin with an (arbitrary) 
orbit segment 70 and its image fi;(7o); the location of 70 is near that of the limit cycle in Fig. [lOl 
The kicked segment clearly folds; indeed, the picture is qualitatively very similar to that of the 
limit cycle case. Note that at Ofb = 1.2, the rotation number of the return map to {O2 = 0} is a 
little below 1, so that <l>t(K(7o)) has an overall, slow drift to the left when viewed in the fixed frame 
[0,1)^. This slow, left- ward drift is not especially relevant in our moving frame (which focuses on 
the movement of $t(K(7o)) relative to that of $^(70)). On successive laps around the torus, the 
orbit in question returns to the part of the torus shown in the figure, and the sequence of actions 
depicted in Fig.[T2]is repeated. We regard this as geometric evidence of shear-induced chaos. 
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Figure 11: Snapshots of the limit cycle and its kicked image in a moving frame. Blue curves: ^^-images of 
7o, the part of the limit cycle between 62 = and 62 = I. Red curves: ^(70) and its images. Purple dot: 
the point on 70 which does not move under k. The parameters are the same as in Fig. [TOl 

We have seen that in the phase-locked regime, the folding of the limit cycle (when the time 
interval between kicks is sufficiently large) can be deduced from the geometry of the strong stable 
foliation. A natural question is: in the quasi-periodic regime, are there geometric clues in the 
unforced dynamics that will tell us whether the system is predisposed to chaotic behavior when 
forced? Since folding occurs in finite time, we believe the answer lies partially in what we call 
finite-time stable manifolds, a picture of which is shown in Fig. [131 We first explain what these 
manifolds are before discussing what they can — and cannot — tell us. 

Fix t > 0. At each z E T^, let V{z) be the most contracted direction of the linear map 
D^t{z) if it is uniquely defined, i.e. if v is a unit tangent vector at z in the direction V{z), then 
\D^tiz)v\ < \D^t{z)u\ for all unit tangent vectors uat z. A smooth curve is called a time-t stable 
manifold if it is tangent to V at all points; these curves together form the time-t stable foliation. In 
general, time-t stable manifolds are not necessarily defined everywhere; they vary with t, and may 
not stabilize as t increases. When "real" {i.e. infinite-time) stable manifolds exist, time-t stable 
manifolds converge to them as t ^ 00. 
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Figure 12: Snapshots of an orbit segment and its image after a single kick in a moving frame, for the system 
(flOl ) in a near-periodic regime. Blue curves: a segment 70 of an orbit and its forward images ^tilo) at 
t = 3.5, 5. Red curves: k(7o) and its forward images. Purple dot: the point on 70 which does not move 
under k. The parameters ai^e ui = 1, V2 = 1-1> Off = 1> and Ofb = 1.2. 

The blue curve in Fig. [13] is an orbit segment of $i. The angles between this segment and the 
time-5 stable manifolds (green curves) reflect the presence of shear. For example, if a kick sends 
points on the blue curve to the right, then within 5 units of time most points on the kicked segment 
will lag behind their counterparts on the original orbit segment — except for the point with 6*1 = 
at the time of the kick. Pinching certain points on an orbit segment while having the rest slide 
back potentially creates a scenario akin to that in Fig. 2; see Sect. II. 2[ One is also likely to find 
shear along the black curve in Fig. [131 a second orbit segment of $4. Whether or not the shear 
here is strong enough to cause the formation of folds in 5 units of time cannot be determined from 
the foliation alone; more detailed information such as contraction rates are needed. What Fig. [13] 
tells us are the mechanism and the shapes of the folds if they do form. Notice also that shearing 
occurs in opposite directions along the blue and black segments. This brings us to a complication 
not present previously: each orbit of spends only a finite amount of time near, say, the blue 
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Figure 13: The time-5 stable foliation of the system (flOl ) in a near-periodic regime. Blue & black curves: 
two orbit fragments. Green curves: time-5 stable foliation. The parameters here are the same as in Fig. [12] 

curve before switching to the region near the black curve, and when it does so, it also switches 
the direction of shear. Finite-time stable foliations for system (flOl) have also been computed for 
t E {3, 5} and a sample of Ofb E (1.1, 1.6) (not shown). They are qualitatively similar to Fig.[T3l 
with most of the leaves running in a northeasterly direction. 

In summary, for t not too large, time-t stable foliations generally do not change quickly with t 
or with system parameters. They are good indicators of shear, but do not tell us if there is enough 
shear for folds to form. For the system defined by (flOl) . given that the finite-time stable manifolds 
are nearly parallel to flowlines and the kick map acts unevenly with respect to this foliation, we 
conclude the presence of shear. Fig. [I2land similar figures for other Ofb (not shown) confirm that 
folding does indeed occur when the system is forced in the near-periodic regime. 

Computation of Lyapunov exponents 

To provide quantitative evidence of shear-induced chaos in the situations discussed above, we com- 
pute Amax- Recall that while periodic kicks followed by long relaxations provide a simple setting 
to visualize folding, it is not expected to give clean results for Amax because of the competition 
between transient and sustained chaos (see Sect. 11.21 ). Continuous-time random forcing, on the 
other hand, produces numerical results that are much easier to interpret. 

Study 4a: Stochastic Forcing. We consider system (flOl) with = I and Ofb E [1.1, 1.6]. The 
forcing is of the form I{t) = a ■ dBt where Bt is standard Brownian motion. 

Study 4b: Periodic kicks. The equation and parameters are as above, and the forcing is given 
hyI{t)=A-Y.nKt-nT). 
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Figure 14: Lyapunov exponent of the system (iTOl ) subjected to white noise forcing. The parameters corre- 
spond to those in Figs. [TO] and [121 respectively. 
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Figure 15: Lyapunov exponent of the system (ITOl) subjected to periodic kicks. The parameters correspond 
to those in Figs. [TO] and [121 respectively. As in Study 1, we show both upper and lower estimates of Amax- 

Summary of Findings. Positive Amax are found for stochastic forcing in the parameter interval 
studied, both before and after the onset of phase-locking at Ofb = a^. For periodic kicks with large 
enough A and T, it appears that A^ax is positive for a fraction of the forcing periods tested, but 
the results are hard to interpret due to the competition between transient and sustained chaos. 

Supporting Numerical Evidence. Fig. [Ml shows some results for stochastic forcing. For Ofb = 
1.47, negative Lyapunov exponents are found for very small amplitudes of forcing, while slightly 
stronger forcing {e.g. a ^ 0.4) is needed before A^ax > can be concluded with confidence. In 
contrast, even fairly small values of forcing seem to lead to A^ax > when an, = 1.2, i.e. in the 
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near-periodic regime. This may be explained by the damping in the limit cycle case, especially for 
larger Ofb. Notice also that in this model large amplitudes of forcing do not lead to larger Amax- 
This is due to the fact that unlike the system in Studies 1-3, a very strong forcing merely presses 
most of the phase space against the circle 9i = 0, which is not very productive from the point of 
view of folding phase space. Fig. [15] shows plots of Amax for periodic kicks. Here, roughly 40% 
of the kick periods T for which Lyapunov exponents were computed yield a positive exponent. 
More generally, we find that Amax > for over 40% of kick intervals T as A varies over the range 
[0.75, 1.5]. See Simulation Details in Study 1. 

Conclusions 

Shear-induced chaos, by which we refer to the phenomenon of an external force interacting with 
the shearing in a system to produce stretches and folds, is found to occur for wide ranges of 
parameters in forced oscillators and quasi-periodic systems. Highlights of our results include: 

(i) For periodically kicked oscillators, positive Lyapunov exponents are observed under quite 
modest impositions on the unforced system and on the relaxation time between kicks (in 
contrast to existing rigorous results). These regimes are, as expected, interspersed with those 
of transient chaos in parameter space. 

(ii) Continuous-time stochastic forcing is shown to be equally effective in producing chaos. The 
qualitative dependence on parameters is similar to that in deterministic forcing. We find that 
suitably directed, degenerate white noise is considerably more effective than isotropic white 
noise (and additive noise will not work). We have also found evidence for an approximate 
scaling law relating Amax to a. A, and a. Other types of random forcing such as Poisson 
kicks are also studied and found to produce chaos. 

(iii) The shear-induced stretching-and-folding mechanism can operate as well in quasi-periodic 
systems as it does in periodic systems, i.e. limit cycles are not a precondition for shear- 
induced chaos. We demonstrate this through a pulse-coupled 2-oscillator system. Chaos is 
induced under both periodic and white noise forcing, and a geometric explanation in terms 
of finite-time stable manifolds is proposed. 

The conclusions in (i) and (ii) above are based on systematic numerical studies of a linear shear 
flow model. As this model captures the essential features of typical oscillators, we expect that our 
conclusions are valid for a wide range of other models. Our numerical results, particularly those 
on stochastic forcing, point clearly to the possibility of a number of (rigorous) theorems. 
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